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Abstract 

The conservation laws for a class of nonlinear equations with vari- 
able coefficients on discrete and noncommutative spaces are derived. 
For discrete models the conserved charges are constructed explicitly. 
The applications of the general method include equations on quantum 
plane, supersymmetric equations for chiral and antichiral supermul- 
tiplets as well as auxiliary equations of integrable models - principal 
chiral model and various cases of nonlinear Toda lattice equations. 



1 Introduction 

Our aim is to present the procedure of derivation of the conservation laws and 
consequently the conserved charges for a certain class of nonlinear equations 
with variable coefficients. In the classical field theory the conserved currents 
and charges follow from the Noether theorem and are connected with the 
symmetry of the action. In the case of linear equations of motion with con- 
stant coefficients the construction of conserved currents by Takahashi and 



Umezawa method pUI can be applied 
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In the previous papers []TU|, [TT], |T^, |TB| we have extended this procedure to the 
hnear equations on discrete and noncommutative spaces including quantum 
Minkowski ||2^, |2^, and braided hnear spaces 0, . It appears 

however that we can consider in a similar way a wide class of equations 
with variable coefficients built within framework of any differential calculus 
with the Leibnitz rule for partial derivatives deformed via the transforma- 
tion operator which is multiplicative and invertible. The range of admissible 
spaces includes the classical space-time with continuous coordinats, the dis- 
crete space, the mixed space with discrete and continuous coordinates, su- 
perspace including space-time and spinor coordinates, quantum Minkowski 
and braided linear spaces (with q- Minkowski as a special case). 
The equations we shall investigate have the variable coefficients which fulfill 
the corresponding restriction (^) including the conjugated derivatives. The 
possibility of derivation the conservation laws for some equations with vari- 
able coefficients on noncommutative spaces was indicated earlier |]12|, but 
the proposed conditions for coefficients were too strong to construct usefuU 
examples. In constrast the new condition (as we show explicitly in the appli- 
cations) appears to be identical with the nonlinear equations of some of the 
integrable models for which we derive the conserved currents via the auxil- 
iary linear equations with variable coefficients. 

Resuming the proposed method can be applied to nonlinear models in two 
ways namely we can consider the equation with nonlinear term free from 
the derivatives (in form of the potential) or alternatively we investigate the 
auxiliary linear equations with variable coefficients for nonlinear integrable 
models. 

The paper is organized as follows: the section 2 contains the description 
of the investigated models and the derivation of the conservation laws. In 
section 3 we extend the procedure to discrete and mixed discrete and contin- 
uous models including both types of derivatives initial d and its conjugation 
. For this class of equations we add the construction of conserved charges 
following the classical method of integrating the time-component of the con- 
served current over the subspace excluding the time-coordinate. This result 
is due to the fact that for the discrete and classical differential calculi the 
definite integral over mentioned subspace is known to commute with time- 
derivative (0). 

The section 4 includes applications of the derived procedure. We start with 
a simple example of an equation of the second order with variable coefficients 
on quantum plane, then we consider the conserved currents and charges for 
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a pair of chiral and antichiral supermultiplts on D=4 N=l superspace con- 
nected via nonlinear equation of motion. This section is closed with a review 
of conservation laws and conserved charges for a range of integrable mod- 
els. We apply our method to the principal chiral model equations and their 
auxiliary linear equations. In this example only the classical commutative 
derivatives are involved; nevertheless it illustrates well the connection be- 
tween the imposed condition @ and the initial nonlinear (in this case the 
principal chiral model) equation. The next model is the nonlinear Toda lat- 
tice equation for which we use our results for mixed discrete and continuous 
spaces in order to write explicitly the conserved currents and charges. 

2 Nonlinear equations with variable coeffi- 
cients and their conservation law 



In the previous papers we discussed the conservation laws for linear equations 



with constant coefficients for discrete differential calculus |jTO|, |TT| as well as 



for a wider category of equations acting on noncommutative spaces namely on 



quantum Minkowski spaces [0, and the braided linear ones fl^, |T^, |T6[ . 
Now we would like to present our results for an extended class of equations 
with variable coefficients of the form: 

A(a)$ = (1) 
A{d)=Ao + ELK...,,d'^'-d''^ (2) 
where coefficients (may be matrices) fulfill the condition: 

"'A,,...,, = (3) 

for / = 1, and with the conjugated derivative described below (PHj). 
We include into considerations the nonlinear equations provided the only 
nonlinear term does not depend on derivatives that means: 

Ao = Ao($) (4) 

As will be shown in the sequel this class of equations contains the equations 
of motion for supersymmetric models in superfield formulation as well as aux- 
iliary linear equations yielded by nonlinear integrable models as for example 
Toda lattice equations or principal chiral model equations (the last being an 
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example built using classical commutative derivatives). 

Our construction presented below in Propositions 2.1 and 2.2 holds for any 

differential calculus with the following deformation of the Leibnitz rule: 

d\f9) = {d'f)9 + {qf)d^9 (5) 
where C is the invertible transformation operator. 

Starting with the classical commutative differential calculus (for which ( is 
the identity operator) we can include the discrete differential calculus with 
derivatives defined by one of the formulas: 

d'fk = (C - l)fk = fk+i - fk (6) 
d^fk = (xk+i - Xk)~^{C - ^)fk = (xk+i - Xk^^ifk+i - fk) (7) 

where k denotes the position on the lattice in the direction i. 
The first one is widely used in the discrete models and measures the difference 
between the values of the function between two points of the lattice in the 
given direction i while the second one (introduced JlOl, |Tl|) is a quotient 
of the change of value of the function over the distance between the two 
neighbouring points of the lattice taken in the direction i. Both types of 
the derivatives obey (|^) with the transformation operator being the shift 
operator along the lattice in the direction i: 

Cfk = fk+1 (8) 

where k in the above formula denotes the value of the i-th coordinate of the 
point on the lattice. 

It is clear that we can also include mixed models depending both on classical 
commutative derivatives with respect to continuous coordinates and discrete 
derivatives (|]) or (^ with respect to lattice dimensions. 
The Leibnitz rule of the form (^ is also characteristic for noncommutative 
spaces. We have checked |l^ that for quantum Minkowski spaces and 
for braided linear spaces (including the q-Minkowski space) it is given by the 
formula (^ with transformation operators determined by their multiplicity 
property: 

Qifg) = iamcfg) (9) 

and the action on monomials of the first order: 

Cjx^ = RfjX'' - {RZ)f (10) 
Qxk = R%xi (11) 
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The first of tlie above formulas is valid on the quantum Minkowski and the 
second on the braided linear space with R being the matrices fulfilling the 
QYBE p4| , p!7[1 . For quantum Minkowski spaces introduced by Podles and 



Woronowicz the R matrix is selfinvertible R^ = 1 and for braided linear 
spaces developed by Majid it is bi-invertible: 

(R-'TX = mR-')^ = m (12) 

Let us notice that the important category of models on noncommutative 
spaces are supersymmetric models in the superspace formulation. In this 
framework the noncommutative space is divided into the classical commuta- 
tive Minkowski space coordinates x (with corresponding indices and metric) 
and the spinor coordinates 6 (we shall not specify below the type of spinors 
- Majorana or Weyl). For all equations within this class the part connected 
with space-time is described using the classical commutative differential cal- 
culus while for spinor coordinates we have anticommuting derivatives. The 
first part of the components of the transformation operator looks as follows: 

= C = C = (14) 

where for the supersymmetric models we have denoted the space-time indices 
as /i, u while we use a and /3 as spinor indices (Weyl or Majorana). The 
only nontrivial part of the transformation operator for the supersymmetric 
models are the components which act between the spinor indices namely for 
the monomials of the first order they are of the form: 

We can extend the supersymmetric transformation operator to an arbitrary 
function by its multiplicity property (|^). 

For all the considered differential calculi (discrete, supersymmetric, quantum 
Minkowski and braided) the transformation operators have one common fea- 
ture namely they are invertible; for discrete models with derivatives (^ or 
d^) the inverse operators are simply backshift operators in the given direction 
while for noncommutative models the inverses are given by their multiplicity 
property: 

C \f9) = (C V)(C '9) (16) 
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together with their action on monomials of the first order: 

(- 'x^ = R]\x'' + Zf (17) 
C" 'xk = R%xi (18) 
where R is the second inverse of the R matrix characteristic for braided hnear 
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spaces 

The inverse operator for supersymmetric models is given by: 

C " = K C " = = (19) 

with the spinor - spinor part defined by the multiplicity (|16D and its action 
on monomialas of the first order: 

g "x, = g -e, = -sp, (20) 

All the described above inverse operators fulfill the condition: 

C^C ' = C7 'a = ^ (21) 

The properties of the transformation operators C and C~ imply the following 
modification of the Leibnitz rule: 

d%Q = i-d^ 'f)9 + fid'g) = / (- * +0"^ 9 (22) 

where the conjugated derivative d'^ is defined as follows: 

at ' ■= -d\- ' (23) 

and we take for the given type of the derivative the connected inverse trans- 
formation operator. We see that after modification we deal with the Leibnitz 
rule where the right-hand side is analogous to the classical differential calcu- 
lus: 



for which the conjugated derivative is given by: d'^ * = —d^. 

The conjugated derivatives form the conjugated equation which for the class 

described by conditions (|l| - H) looks as follows: 



-t, , , X ^ Ml Mi 



Aid ) = Ao(<f ) + -d {g "^-C, "^A.,....J (24) 



=1 



The following propositions hold for equation (|1| - 0) and are the extension of 



results derived in [11, 12 
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Proposition 2.1 The unique solution of the operator equation: 

" +Q^) o r,{d, d^) = A{d) - A(d^) (25) 



in the class of polynomials of derivatives d and d for the equation operator 
A fulfilling is of the form: 

^,{d, 5^) = (C "A.) + (26) 

1=1 k=0 

Proof: 

Let us explain that the "o" operation describes the way the noncommuting 
derivatives work on monomials of derivatives d and d, namely: 

-t M 



(- d +0'') o [ui,...,ui]a{x)[pi,...,pk] :-- 



-[ui, fi]a{x)[pi, pk] + [iyi,...,ui]d^'a{x)[pi,...,pk] (27) 
where the following notation for monomials was used: 

[pi,...,Pk]:=dPK..d'"^ (28) 

-t 1/1 Ul 



[iyu...,ui]:=d ...d (29) 

The further calculations base on the assumption that we consider the general 
operator of the order — 1 with respect to the derivatives and on the asso- 
ciativity of the algebra of derivatives. The explicit solution of the condition 
( [25| ) is enclosed in Appendix A. 

We modify the F operator due to the deformation of the Leibnitz rule (^ ^2D 
and obtain the operator V: 

f,(9,aVc^ '(C-"A.)+ (30) 

1=1 k=0 



which for a pair of arbitrary functions / and g is connected with the F 
operator by the equahty: 



(31) 



The immediate consequence of the Proposition 2.1 is therefore the construc- 
tion of conserved currents yielded by the following statement. 

Proposition 2.2 Let us assume that function $ is an arbitrary solution of 
equation ([|, ^ with coefficients fulfilling that means: 

A(a)$ = (32) 

and function $' solves the conjugated equation with the operator of the form 

$'A(a ) = (33) 

Then 

j^ = $'f^(a,aV (34) 

where the operator f ^ is defined by (^^, is a current which obeys the con- 
servation law: 

J29''J^ = (35) 

Proof: 

This is a corollary from the (0) property of the T operator, from the operator 
equation (^) for F and finally from the fact that functions $ and $' fulfill 
the respective equations ([l|, |33|) . 

Let us notice that the auxiliary conjugated equation for nonlinear model is 
now the linear one with respect to the solution 

Ao($) + 5^ a ... a C,- "'-C, "^A«....«J = (36) 

The solution $' depends on the solution of the initial equation $ which 
defines the potential for the conjugated equation Ao($). 
Now the interesting question arises whether the presented construction which 
works on a wide class of spaces can be extended to equations with operators 
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including both d and 5^ derivatives. As we know from discrete models (see for 
example P, the equations depend explicitly on the forth- and backshifts 
along the lattice. Also the equations of motion built within the framework of 
generalized difference derivatives (which was discussed in ||^, |lT| ) include 
both initial and conjugate derivatives as a consequence of the minimal action 
principle. The method of derivation of the conserved currents for the discrete 
and mixed models shall be described in the next section. It is based on the 
fact that acting on functions of commutative coordinates the operators (^ 
and d^) yield two symmetric formulas for the Leibnitz rules (both for initial 
and conjugated derivative). 

3 Conserved currents and charges for discrete 
and mixed models depending on the d and 
5^ derivatives 

Our aim is now to extend the construction described in Propositions 2.1 and 
2.2 to the equations depending on the initial and conjugated derivatives of 
the form: 

[A(a) +A(<9t)]$ = (37) 
A(a)=Ao + E;iiA,,...,,a^-...a^' (38) 

A(at) = Ao + e£iAmi...m.5^ ^'■■■d'^ (39) 

Following the previous considerations we restrict the class of equations to 
variable coefficients fulfilling the extended version of the condition (^: 

at '^^A,,...,, = d^'^K.-,, = (40) 
for all / = 1, and A; = 1, iV. 

We include the nonlinear equations provided the nonlinear terms do not 
depend on the derivatives that means they have a form of potential: 

Ao = Ao($) Ao = Ao($) (41) 

As we shall see in the subsequent sections the class of equations fufiUing (^Of ) 
contains the widely studied in different context linear auxiliary equations 
connected with nonlinear integrable models such as principal chiral model 
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and equations of nonlinear Toda lattice models. 

The important factor in our calculations is the fact that for discrete (or 
mixed - classical and discrete) differential calculus we can write the modified 
Leibnitz rule symmetricaly in the following form: 

5"(C ^f)9 = i-d^ ^f)9 + fd^g = fi-d^ ^ + d^)9 (42) 
"(Cf /)^? = {-d^f)9 + fd^ ^9 = fi-d^ + ^)9 (43) 

The second of the above equations indicates that the statement analogous 
to the Proposition 2.1 holds for the operator A((?^); namely we shall prove 
that the operator equation given below has the unique solution by virtue of 
Proposition 3.1: 

J](- 7 n ° f ^(5^ 9) = A(9t) - A(a) (44) 
where the conjugated operator A(5) looks as follows: 

N 

A(a) = Ao($) + J]9 ...a (C-CA,,....,) (45) 
1=1 



Proposition 3.1 The unique solution of in the class of polynomials 

of derivatives d and d'^ for the equation operator A with coefficients 
fulfilling 1 ^4^ ) is of the form: 

r,{d\d)= (46) 

N-l I 

(CAa) + E E - (CC-CA,,...,,.,,^,..,^) at ^^^K..d^ 

1=1 k=0 

Proof: 

The formula P6| ) is the immediate consequence of the proof of Proposition 
2.1 enclosed in Appendix A as the following connections between derivatives 
and transformation operators hold: 

(C^)'^ = C (48) 

(at)t 7 = -{-d^g °)c = d-y (49) 
10 



We notice that C~ is the transformation operator for the derivative d\ the 
operator C being its inverse. Therefore all the calculations from the proof 
of the Proposition 2.1 can be repeated with suitable replacements. In this 



way we obtain the formula for the T operator p6|) which is simply the (p6D 

operator with new derivatives 9^ and d and the operator C acting now as the 
inverse transformation operator. 

We modify the operator F similarly as in (RDI) and obtain the operator F: 



^;.(9^a)=c^(C^Aa)+ (50) 



N-1 I 



1=1 k=0 

which in turn fufills the equality: 



for a pair of arbitrary functions / and g. 

We can use the constructed operators F and F to derive the conserved cur- 
rents by the following proposition: 

Proposition 3.2 Let us assume that the function $ is an arbitrary solution 
of equation with coefficients fulfilling (^J\ ) which means: 

[A(a) +A(9t)]<l> = (52) 



and the function $' solves the conjugated equation with operators (\2^, \43i ) 



^'[A(d^) + A(d)] = (53) 

Then: 

J, = $'f J, = $'r^$ (54) 



where the operators F and F are given by ( y^JJq j is a current that obeys the 
conservation law: 

5^a^J^ + 5^at ^^J^ = (55) 
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Proof: 

The conservation law (|55D is implied by the properties (|31], 0) of the applied 
operators F and F: 

d^^'t^^ + /^^T^* = (56) 

( $^(- oF^(a,aV$^(- d' +at ^) of^(at, a) J $ = 

\ fj- fj- / 

$' (^A(a^) - A(^) + A{d) - A(d^)^ $ = 

Corollary 3.3 T/ie components of the conserved current for ( [5^ fulfilling 
the condition can &e written also in the form: 



j; = <l>'F^<l> - C^- ^ ($'F,$j (57) 
They obey the conservation law including only the d derivatives: 

J29'J', = (58) 

provided $ is the solution of and the function $' solves the conjugate 
equation with operators ('2^,Jd{). 

Let us notice that the presented construction holds when we work within 
the framework of classical commutative differential calculus (for which the 
conjugation means = —d and the transformation operator is the identity) 
as well as for mixed models where part of the coordinates is continuous and 
the other part is discrete with suitable difference derivatives in the form (^ 
or (0). 

In the sequel we shall use the construction (^, ^ ^3|) in order to derive 



the conserved currents and in consequence the conserved charges for mixed 
discrete nonlinear Toda model and for double discrete Toda model as well. 
As we know from the classical field theory we can use the conserved currents 
in construction of the conserved charges. This is also the case for our dis- 
crete and mixed discrete models with derivatives defined by (^ or (|^). We 
shall integrate the time component of the conserved current over subspace 
excluding the time-coordinate. We use in derivation the conserved currents 
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given by Proposition 3.2 or accordingly by the Corollary 3.3 if the conjugated 
derivative with respect to time-coordinate appears in the model. Integrating 
we must remember to take the corresponding integral namely for continu- 
ous coordinates we understand J dxi as the Lebesque integral while for the 
discrete derivative (||) the definite integral is given by: 

/+00 
dx, := J2 (CO' (59) 

A;=— oo 

If the equation includes the discrete derivative of the type (|^ we use the 
definite integral in the form: 



-co 



dxr.= ^ (4+1 - 4) (0)' (60) 



In the above formulas for discrete integrals the Q is the shift operator in the 
direction i along the lattice. 

Let us denote the integral over subspace with excluded time-coordinate as 
j^^y Then the conserved charges for discrete and mixed discrete models with 
variable coefficients can be derived using the following proposition: 

Proposition 3.4 Let us assume that in the model described by the equation 



the conjugated discrete derivative with respect to the time- coordinate 



does not appear. Then the charge 

Q = I Jt (61) 



sub 



where Jt is the time component of the conserved current described in the 
Proposition 3.2 is conserved: 

d'Q = (62) 

Proof: 

The conservation of the charge ( |6TD is implied by the conservation law from 
Proposition 3.2 namely: 

d'Q= [ d'Jt= I -J2{d'Jk + d^ 'Jk)= (63) 

J sub J sub j^^^ 

= boundary terms = 
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provided the currents vanish at the infinity in the corresponding space-hke 
dimensions. We have used here also the property of the discrete and mixed 
discrete differential calculus: 



a* / = / a* (64) 

J sub J sub 



Let us point out that the equality ( |51D does not apply in general to non- 
commutative spaces (the superspace being an exception) and this fact is the 
major difficulty in construction of the conserved charges for an arbitrary 



model on noncommutative space |T^, 116 1 



If in the equation (pTQ both initial and conjugate time-derivatives appear 
we derive the conserved charge using the Corollary 3.3 namely the following 
statement is then valid (with proof being the copy of the above calculations 
for Proposition 3.4): 

Corollary 3.5 Let us assume that in the model described by the equation 
( |5^ both discrete derivatives 5* and d'^ * appear. Then the charge 

Q = I J't (65) 

J sub 

where J[ is the time component of the conserved current described in the 
Corollary 3.3 is conserved: 

d'Q = (66) 



4 Applications 

4.1 The second order equation on quantum plane with 
variable coefficients 

Let us start with the simple example of the second order equation on quantum 

plane with coefficients depending on coordinates x and y. 

The commutation relations for derivatives and coordinates look as follows 



0: 



yx = qxy d^ff" = q ^d^^d^ (67) 



14 



The Leibnitz rule 
matrix: 

R 



in differential calculus is defined by the following R 

R! = q-^R 

















<1 













q 

















Form the above matrix we can deduce the explicit action of the transfor- 
mation operator C, on the monomials of the first order (|TT]). It yields the 
following formulas for the inverse transformation C~: 

C = q~^^ C ""y = q'^y Cy = iq"^ - 1)?/ C ''y = o (69) 



Q'y 



q ^y 



q yy = o 



'X 



(70) 



The conjugated derivatives are defined using the inverse transformation op- 
erator: 



at 



Q] y = _QX^- y 



(71) 



Now let us discuss an arbitrary equation of the second order with the coef- 
ficients depending on x and y and check which of them fufill the condition 

i): 

A(9", d^) = A,,a"a" + Aya^dydf + Kydfd^ + A^^a^a^ (72) 
The condition (^ for our simple case reads as follows: 



"A,, + at yky, = 



at ykyy + at "a 



xy 



Let us assume A 



xy 



^yx 



0. Then we can show that A^ 



(73) 
= (j){y) and 



Kyy = ifj^x) obey the conditions (|), namely: 



at ^•A.,. = (-a^'c ^-a^C" n<P{y) = -^^'C Xi/) = -d'<P{q^'y) = o (74) 
at yAyy = (-a^c ^-a^g ^)^(a;) = -o^c; "^H^) = -oy^j^q-^x) = o (75) 

Now we can apply our construction to obtain the F operator: 



-t X 



a (l^{q-^y) + ct>{q~'y)d^ 



-t y 



a i}{q-^x) +i>{q'^x)dy (76) 



It is easy to check that this operator fulfills the operator equation (P^): 

{-d^ % d^) oT^ + {-!) % dy) oTy 
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where the operator A{d) gives the conjugated equation: 

A(a) = d d + a a ^piq-''x) (78) 

The modified f operator is given by the formula: 

fx = ^<P{q~'y) + Cx. (79) 

+ "^y %(g-'x) + Cj, V(g"'2:)9?' (80) 

Having the exphcit form of the F and F operators we can construct the 
conserved currents: 

J,. = $T,<I> = $T,$ (81) 

where the functions $' and $ solve the respective equations: 

(^a^ "^(g-'y) + a^ ""^(g-^s)^ = (82) 

(0(?/)a^a^ + ^{x)dydy) $ = (83) 

According to the Proposition 2.2 the above current obeys the conservation 
law: 

+ dyjy = Q (84) 

provided functions $' and $ are the solutions of the corresponding equations 

4.2 Nonlinear equation of motion for chiral and an- 
tichiral supermultiplets 

The supersymmetric models in the superfield formulation yield interesting 
examples of the equations of motion with coefficients depending on variables 
H, p5[ . Let us recall that in such a framework the fields depend on superspace 
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coordinates including the space-time and spinor variables while in construc- 
tion of the action the covariant derivatives are used which explicitly depend 
on spinor coordinates. 

We shall study in this section the D=4 N=l chiral and antichiral superfields 
obeying the nonlinear equation of motion resulting from the supersymmetric 
action: 

1 = 1 d^xdpesm^ + j d^xcfe (^$2 + + j d'xd^e (^$2 + 1<|3 

(85) 

The integration over spinor variables 6 and 6 can be expressed in terms of 
the covariant derivatives: 

Jd^e = -lD^ Jd^9=-lD^ (86) 
/ dWe = ^Om^ = ^D^D^ = ^D'^D^D^ = ^D^D^'^ (87) 

The covariant derivatives used in supersymmetric models are built from the 
basic derivatives with respect to the space-time and spinor coordinates: 

D^ = d^ + la^^Pd^, = e^fDf, (88) 

D. = _a. _ zO^a'^^d^ D° = e^^D^ (89) 

Due to the fact that chiral $ and antichiral $ superfields fulfill the following 
condition: 

= D^^ = (90) 

we obtain from the action (|85D the equations of motion in the form: 

= (91) 



jD"^ m + \ f $ 



■m + \D'^ ) \ ^ 

We have written the equations for the superfields $ and $ in the matrix 
form. In this example we deal with the coefficients of the equation depending 
explicitly on 6 and 9 variables due to the form of the covariant derivatives 



= D"D^ = -d'^da + 2ie'^a'^-df,d" - 9^0 (92) 
= DaD^ = -4a" - 2i9'^a''-,dJ^ + 9'^U (93) 
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The second important feature of the chiral-antichiral superfield equation of 
motion is its nonhnearity for (7 7^ 0. However the nonhnear term does not 
depend on the derivatives so we can apply our method of derivation the con- 
servation law provided the condition (|^) for the kinetic part of the operator 
is fulfilled. Let us extract the ( and C~ operators from the Leibnitz rule for 
basic derivatives: 

d^{fg) = {d^f)g + fd^^9 (94) 
d'^{fg) = id^f)g+iQf)d^g (95) 

d''{fg) = {d''f)g+iCif)d^g (96) 

It is clear that Cj^ = 6^ and the following components of the transformation 
operator vanish: 

C^ = C = o C = C = o C = C = o (97) 

The remaining components of the ( operator are defined by its action on the 
monomials of the first order in superspace variables: 

C^x, = 5|a;, g9^ = -6''^e^ ^9^ = -Sp^ (98) 

The analogous formulas hold also for dotted indices. Hence we get for the 
inverse transformation operator: 

C/3 = SpXi, 6^ = —5^9^ C,p 9^ = —5p9j (99) 

with the corresponding expressions with dotted indices of the same form. 
Now having the explicitly derived inverse transformation operator we arrive 
at the conjugated derivatives given by the formulas: 

5t ^ = -d^" " = -df^g " ^ = -d^q " (100) 

Using the conjugated derivatives we obtain the condition (H) written for our 
equation for the chiral and antichiral superfields: 

at °A,^ + at '^A^^ + at -a,^ = o (loi) 
at -A,^ + at '^A^^ = (102) 
at -A + at '^A . = (103) 
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One can easily check that the above conditions are fulfilled using the explicit 
form of the coefficients from equation (|9T|): 

Kf^ = = (104) 
Aa/3 = -e„^ ( ) 

= -^aB ( n ^i) (106) 



=^a,= (^ q"""^ ° ) gu, (107) 

A^id = Aq^ ~ —-d^a'^p- ) ^'^'^ (lOS) 

K^=y I _i_02 ) 9,u (109) 

The next step is therefore the application of the general formula (^) for the 
r operator to our model: 

= 'a.m + Am.<9" (110) 

r„ = -d^ + + A„^a^ - K^.^^' (111) 

= -a^ + "a^^ + A^^a^ - A^^a^ (112) 

The operator F differs from the F by insertion of the C~ operator in the 
middle of the monomials of derivatives: 

fM = r^ (113) 

tp = -^^^Xfi "a^„ + a^ "a^„ + C/3 ^A^^a^ - C/3 "a.^S'^ (114) 

f ^ = -a^ C/3 "a^^ + a^ "a^^ + C/3 ^A^^a^ - C/j ^A^^a^^ (115) 



We apply the obtained operator F to derive the currents: 

$ 



•^'O ( I ) (116) 
J„ =($',$') fa ( I ) (117) 
J, = ($', $0 f d f I ) (118) 
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which obey the conservation law: 



(119) 



provided the fields used in construction fulfill the corresponding equation of 
motion (|9TD and the pair of superfields $' and $' its conjugated version: 



Let us notice that this equation is linear with respect to superfields $'and $' 
and its solution depends on the potential given by the chiral and antichiral 
superfields $ and $. 

Due to the properties of covariant derivatives built form the self conjugated 
operators we conclude that the conjugated set of equations for superfields 
($', $') has at least following solutions: 

• for the case (7 7^ we can take $' = $ and $' = $ where $ and $ are 
solutions for (|9l|), 

• when g = Q the superfield is an arbitrary solution of chiral and 
of antichiral part of the set of equations (^Tj). 

The above solutions give in turn the explicit form of the conserved currents 
( |116| - |118| ) where we replace the multiplet (i>', $') with ($, i>) ioi g ^ Q and 
with (5$, 5$) for g = Let us notice that for linear model 51 = we 
obtain the full set of conserved currents connected with symmetry operators 
of chiral-antichiral supermultiplet equation: 




(120) 




(121) 



(122) 



(123) 



where the symmetries 5 include supersymmetric transformations: 



a 



(124) 
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and the operators from Poincare algebra: momenta, angular momentum and 
boosts for four-dimensional Minkowski space. 

We can use the obtained conserved current ( |116| - |118| ) to construct the 



integral of motion by integrating the time-component of the current ( |116|) 
over the subspace: 

Q = J £x(fe(feJo (125) 
as the following equality holds: 

= J £xd^ 9(1^9 d^Jo = (126) 

j d^xd'^9d^9 {-d^Jk - a° Ja - d^J^) = boundary terms = 

Let us point out that the developed method allows immediate construction 
of conserved charges. The obtained charges are supermultiplets built from 
component charges which are also conserved separately. The application 
of our method to the equations obeying the condition (|^) is an alternative 
to the procedure used in supersymmetric models (see for example [Q for 
chiral superfields and ^ for supersymmetric principal chiral model) where 
the conservation laws include covariant spinor derivatives D and D instead of 
basic derivatives d'^, 9" and 9" which we have used. The consequence of the 
fact that the time-derivative does not appear explicitly in the conservation 
law is an additional procedure required to derive the component conservation 
laws and conserved charges. 



4.3 Nonlinear integrable models and their consistency 
condition 

4.3.1 Principal chiral model equation 

We shall discuss the application of the construction given by Propositions 2.1 
and 2.2 to the classical problem of nonlocal conserved currents and charges 
for the principal chiral model. The model is built using the commutative 
differential calculus, nevertheless its linearized version yields an interesting 
example of the linear equation with coefficients fulfilling the condition (^. 
Dimakis and Miiller-Hoissen p, ^ treated this case within their framework of 
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the gauged bi-covariant noncommutative differential calculus and obtained 
the linearized version of the equations in the form: 



(127) 
(128) 



where x is an x iV matrix of smooth functions depending on t and x while 
the field g denotes the principal chiral field and the initial nonlinear equation 
of the principal chiral model looks as follows: 



d^ig-^d^'g) = d\g~'d'g) 



(129) 



In the paper the integrability condition for the set ( [127|J128|) was used and 
it appeared that it is first order in derivatives therefore yields the conservation 
law with currents linear with respect to the field x ^^^^ depending on the A 
parameter. 

We propose in addition to consider the consistency condition derived from 
the set ( |12?| , |128| ) which has the following form: 



[d^ff" - d'd' + g-\d''g)&' - g-\d'g)d'] x = 



(130) 



As we see the operator of this equation is free of the parameter A and the con- 
dition becomes identical with the principal chiral model equation ( |129| ): 

"A, + 'At = d'g-^d'g) - d'g~\d'g) = (131) 

the remaining coefficients A of the operator A{d) being constant. 

The conjugated equation for ( |130|) together with connected conjugated set 

for ([1271 , [1281) look as follows: 



X' 



d d — d d 



- d g'\d^g)+ d g-\d'g) 
X[-d-x' + ^g-\d-g)] 







-d-^ = \[-d'x' + x'9-\d'g)] 



(132) 

(133) 
(134) 



We work in this example with classical commutative basic derivatives so the 
transformation operator is simply the identity operator and the operators F 
and r coincide: 



tx = T, = g-\d-g)- d 
ft = Tt = -g-\d'g)+ t 



(135) 
(136) 
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Our consistency condition allows the construction of the conserved current 
according to the Proposition 2.2. The conserved current is given by the 
derived operator F and the solution x of the equation ( |130| ) while x' denotes 
the solution of the conjugated equation ( p.32| - |134| ): 



J, = x'r.^x Jt = }cTt5x (137) 



where we have denoted as 5 the symmetry operator of the equation (|130|) . 
The operator 5 can be derived using the extended vector fields method de- 



scribed in |21|. We have checked that the following operators transform 



solution X of consistency condition ( |130|) into solution 5x'- 

5' = C^ 5' = X^ S' = 9^'S- S'={d^)-\9-'d'9)-^ (138) 
dx dx dx dx 

where C is a constant matrix N x N. 

The above symmetry operators of the equation ( p. 3 Op imply the following 
currents: 

j; = x'r,c Jl = x'r,x 4 = x'r,g-' 4 = x'r.id^'ig-'d'g) 

(139) 

which are conserved: 

ff'Ji + d'Ji = (140) 

thereby yielding the conserved charges after the integration of the time- 
components of the curents (|139|) with respect to the space variable x: 

Q' = Jdx x'^tC (141) 

Q' = Jdxx''^tX (142) 

Q' = Jdx x'Ttg-' (143) 

= J dx x''^t{d'^y\g-^d'g) (144) 



where the operator Tt is given by formula (|136D . 

Now following p| we can expand the fields x ^"^^ x' terms of powers of 
the parameter A and obtain the infinite set of conserved charges connected 



with the nonlinear equation ( |129| ) of the principal chiral model: 

oo oo 



m=0 m=0 
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and (r5^,134|) we arrive at the following set of 



Using equations 

equations for coefficients x*'™'' and ;^'(™) (we assume after 0] that x 
x'^^^ and x'^^^ a-^e the unital N x N matrices): 

d'x^^^ = [d^ + g-\d^g)]x^"'-'^ 
which give after solution: 



m > 2 

m > 2 

m > 2 

m > 2 



(2) 



(9-)-i^-i(9*<7) 



'X 



'(2) 



X 



(m) 



[(9")-i(a* + (7-^5*^7)]"" (s^O^'^^^H^*^?) 

m-2 



(0) ^(1) 

(146) 
(147) 
(148) 
(149) 



(150) 
(151) 

(152) 



Inserting the solution ( |151| , |152D into the conserved currents ( |139| ) we obtain 
after expansion with respect to the parameter A the infinite tower of con- 
served charges for the field g of the principal chiral model: 



go (m) = J^^ x'^'")rtC 
Q2 (m) = J(l^ x'^^^r^C/-! 

Q3 M = x'^™)rt(a^)-i^-i(a*^) 

4.3.2 Generalization of the principal chiral model equation 



(153) 
(154) 
(155) 
(156) 



We can easily extend the results of the previous section to the case of the 
principal chiral model with the fields depending smoothly on variables t, x 
and y. Let us recall the linearized version for N x N matrices x with entries 
depending on the mentioned variables H: 



d^X = A(5^ + ^ \d''g))x 
dyx = X{& + g~\d'g))x 

with the conjugated set looking as follows: 

d^x' = A [d^'x' - x'g'^id^'g)] 

dyx' = X[d'x'-x'9-\d'9)] 



(157) 
(158) 



(159) 
(160) 
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The integrability condition of the set ( |157| , p.58| ) yields the currents discussed 
in 0. Following the previous example we propose to investigate the consis- 
tency condition for ( |157| , |158D which becomes: 



^d^gy _ a* a* + g-\d^g)dy - g-\d'g)d'] x = 
with its conjugation in the sense (^) of the form: 



^d-gy - d'd'] x' + x' 



(161) 



(162) 



The principal chiral model equation coincides with the condition (|^) for vari- 
able coefficients: 



+ 'At = dyg-\d''g) - d'g-\d'g) = 



(163) 



We can write down the components of the F operator taking into account 
that due for the fact of the basic derivatives being commutative, the trans- 
formation operator is the identity so the operators F and F coincide: 



d +0^ 



Ty = Ty=g-\d^g) 



g-Hd'g)+ d -a* 



(164) 
(165) 
(166) 



The current built using the above operators and the solutions x ^-iid x' of 
the respective equations (|157| - |16CI| ) looks as follows: 



4 = x'rjx Jy = xry6x Jf = x'^tSx 

with S being the symmetry operator from the set: 

S' = C^ S' = x-^ S' = 9-'^ S' = idy)-\g-'d'g)^ 
ox ox ox ox 



(167) 



(168) 



where C is a constant N x N matrix. 

According to the Proposition 2.2 the currents (|167|) obey the conservation 
law: 

(169) 



d'Jl + d'Jl + d^Ty = 
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Thus the currents ( p.67| ) produce the infinite tower of conserved charges sim- 
ilarly to the previous case. After expansion of the fields x ^i^nd in terms of 
powers of the parameter A ( |145| ) we obtain: 

go M _ J dxdyx'^"''^TtC (170) 
Qi M _ j^^^y ^™ ^ x'^'^)!^^'"-') (171) 

Q2 (m) ^ J ^^^y ^'{m)Y^g-l (172) 
q3 (m) _ J ^^^^ x'^™)rf(a^)-l^-l(a*C/) (173) 

where the explicit expressions for the component fields given by the set of 
equations ( p.57| - |16CI|) look as follows ( with x^^\ X^^\ x'^^^ x'^^^ being 
the unital N x N matrices as before): 

^(2) ^ _^,(2) ^ (a.)-i^-i(a*<?) (174) 

^ + g'\d'g))r~^ (dy)-^ g-\d' g) m > 2 (175) 



^/(m) ^ (•_l)m-2^/(2) 



-9 +(7-i(a*(7)') (a^)-^ 



m-2 



m > 2 (176) 



4.3.3 Nonlinear Toda lattice equation 



The interesting case of application of the Proposition 3.2 is the linear set 
of equations which originates from nonlinear Toda lattice equation. The 
set of auxiliary equations derived in the paper from gauged bicovariant 
differential calculus reads as follows: 

Xk = A(e^^-i-''^Xfc-i - Xk)C (177) 

Xk+i -Xk = -Kxk + qkXk)C (178) 

where we have used our notation for the transformation operator which for 
the mixed model like the one considered in this section is simply shifting 
along the lattice C/fc = fk-i with respect to the space-like dimension. Again 
we shall discuss the consistency condition which for the set of equations 
( |177| , |178D becomes an equation of second order including classical derivative 
with respect to the time-variable and the discrete one with respect to the 
space-lattice variable: 

Xk + qkXk - e't'-'-'KXk-i - Xk) - {Xk+i - Xk) = ^ (179) 
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Let us notice that this equation is free of the parameter A and the operator 

C- 

We can rewrite this equation using the notation: 



Xk-i - Xfe = (C - l)Xfc = dxk Xk+i - Xk 

Qk-i -qk = iC- i)gfc = dqk f = a*/ 



(C- - l)Xk = d^Xk (180) 
(181) 



and it has the following form: 



182) 



In this form we see clearly that the equation involves both d and derivatives 
so we should follow the construction for discrete models given in Proposition 
3.1 and 3.2. 

Let us observe that the discussed equation fufills the restriction (^ and it 
coincides in this case with the nonlinear Toda lattice equation: 



Qk-l—Qk 







;i83) 



We construct the F operator according to (|2^ , ^ |50D taking into account 
that the transformation operator acts as follows: 



Ctfkit) = fk{t) 



Qfkit) = fk-i{t) Cfkit) = fk+i{t) (184) 

(185) 



Tt = -d +d' + id'q) 



(186) 
(187) 



The components of the f operator yield the corresponding components of 
the current: 



Jf = 6x'TtX 4 = Sx'T.X Ji = <^x'f .X 



;i88) 



which obey the conservation law including the conjugated derivative: 

aVf + dJi + d^Ji = (189) 



provided the field x fufills the consistency condition ( |182| ), the field x' its 
conjugation: 

X' tt -t {d'q)-l) e'^^'^'^-d =0 (190) 
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and the operator S is now the symmetry operator of the conjugated equation 

m- 

We proceed taking integrals with respect to the discrete spatial variable (we 
have denoted it ) and obtain the charges: 

Q^^ji.J.^j,.,Ax (191) 
which are also conserved: 

a*Q = (192) 

In the above formula we understand the discrete definite integral j dx as 
given by (|59|) . 

Analogously to the examples for principal chiral model we can check that 
the following operators transform solution of the equation ( |19(]| ) into the 
solution 5x': 

5° = c^ 5^ = 6"^ 5^=[-d'\d'q)+t]^^ (193) 
dx' dx' dx' 

where the indefinite discrete integral is given by: 

^ oo 

5-^ = T3I = -E^' (194) 



k=0 

and c is a constant. 

The derived symmetries of the equation ( |182| ) lead to the following expres- 
sions for charges: 

QO = j dx cTtX (195) 

= j dx e'iVtX (196) 

Q^ = Jdx [-d-\d'q)+t]rtx (197) 
After expanding the solutions x with respect to powers of the operators AC 

and x' ill terms of A C we shall be able to rewrite the expressions ( |191| ) as 
infinite set of conserved charges: 

go(m) = Jdx cFtx^™^ (198) 

gi(m) = Jdx e^r^x^™^ (199) 

Q2{m) = Jdx [-d-\d'q) +t]TtX^"''^ (200) 
where the components x*-*"-* read as follows: 

= [{d^)-\d' + (d'q))]""'^ {d^)~\d'q) m>2 (201) 
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4.3.4 Nonlinear Toda lattice equation - generalization 

Our next example is the generalization of nonlinear Toda model sketched 
in The functions depend now smoothly on one time variable t, one 
spacelike variable y and one discrete space-like coordinate x. The bicovariant 
differential calculus is in this case defined as follows 0: 



where X and Y are matrices with entries being smooth functions of t and y 
and discrete in the coordinate x. The condition 6 A = leads to the following 
generalizd Toda model nonlinear equations: 



Sf=[S-\f]r-fC 



df = fr+[S,m (202) 
(203) 



A = Xr+ {Y ~1)SC 



X'k — Yk — Yj. 



-1 



Yk — YkXk 



XkYk 



(204) 




and looks as follows: 




(206) 
(207) 



We can rewrite the above set of equations using our notations: 



5 = C 

Cxk{t,y) = Xk+i{t,y) 




(208) 
(209) 
(210) 
(211) 



C Xk{t,y) 



Xk-i{t,y) 





X = d^X 



X = d\ 



and they are given by the formulas: 



dyx = \[dx + {y -mx)]C 



(212) 
(213) 



The consistency equation for the above set of equations reads as: 



[9^9* + Xdy -Yd-d^]x = Q 



(214) 
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Let us check the condition (j^) for the variable coefficients of the operator of 
the equation: 

^Ay + d^A^ = -d^X -Y + CY = (215) 

As before our equation fulfills the restriction due to the fact that X and 
Y obey the generalization of the nonlinear Toda model ( p04| ). The condition 
(|40| ) also holds as the coefficient for the derivative is constant. 
Similarly to the previous examples we can apply the Propositions 3.1 and 3.2 
for the discrete models depending on initial and conjugate discrete deriva- 
tives. In this way we obtain the T operator in the form: 

rt = r, = + dy) f, = r, = + a*) + x (216) 

= -C iCY) L = -C (217) 

We use the above components of the T and T operators to construct the 
respective components of the currents: 

= Sx'tx Jy = Sx'TyX (218) 

4 = Sx'T.x Ji = Sx'tx (219) 

where we have denoted 6 as the symmetry operator of the conjugated equa- 
tion ( |220| ) and x is the solution of ( |212| - pi4| ) while x' solves its conjugation: 

[-a^a* + d]x' + idyx')x + id^x')iC-y) = o (220) 

As the currents fulffil the conservation law: 

d'J^ + d^fy + dJi + d^Ji = (221) 
they yield the conserved charges (with the integral J dx in the sense of (|59|)): 

= j dydx Jf = j dydx 6x'ttX (222) 

Following the procedure applied earlier we expand the solutions x ^-nd x' 
terms of powers of the operator AC and AC respectively: 

oo oo _ 

X = Y.^Kr^'-'^^ x' = E(^?~r^'^™^ (223) 

m=0 m=0 
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After assuming the first components in the form of the unital N x N matrix: 

^(0) = y(0) ^ ^(1) ^ yd) ^ 1^^^ (224) 

we obtain the following equations for the subsequent components: 

g]^(-m) ^ ^gt ^ m> 2 (225) 

= [5 + (y _ l)(];^(™-i) m > 2 (226) 

Inserting the expansion (|223|) into the expressions for currents and charges 
we arrive at the infinite towers of the conserved currents and charges: 

4 ^'"^ = ^x'r.x^") (227) 

q5 (m) ^ J ^^^y ^yftx^""^ (228) 
with the following explicit expressions for the components of the fields x'- 

= = (at)-ix (229) 

_ [(5t)-i(5t + X)]™-2(5t)-ix m > 2 (230) 



and the initial solutions for the conjugated equation ( |220|) related to the 
symmetry operators in the form: 

= C 6'x' = d-'X - y (231) 

with C being a constant N x N matrix. 



4.3.5 Double discrete nonlinear Toda lattice equation 

Let us end the review of the integrable models with fully discrete Toda lattice 
equation. The set of auxiliary linear equations given in has the following 
form: 

Xk{n) - Xk{n - 1) = -^[9k^{n)gk+i{n)xk+i{n) - Xk{n)] (232) 
Xkin) - Xk-iin) = --fc[gi:\n)gk{n + l)xkin + 1) - Xkin)] (233) 

We rewrite the above equations using our notation: 

fk{n) - fk-i{n) = -d^fkin) Q fk{n) = h+i{n) (234) 
fk{n) - Un - 1) = -d'Un) Q fk{n) = /,(n + 1) (235) 
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Then the set of equations ( |232| , |233| ) looks as follows: 

-d\ = -l[9-\C^9mx) - X] (236) 
-d-x = -lc[g-\Q9){Cx) - X] (237) 

The consistency equation for the above set reads as: 

- 9-\Q9)d'^ 'X - d^X + c^9~\Ct9)d^ \ + = (238) 

The condition ( ^O]) for our example is fulfilled: 

^A^ + at 'At = (239) 
d^A^ + d'At = -d^g'HC9) + c^d'g-^Cig) = (240) 

provided the double discrete nonlinear Toda lattice equation is valid for the 
field g : 

d''9-\C9) = c'd'g-\Qg) (241) 

We construct the conserved currents following Propositions 3.1 and 3.2. As 
both initial and conjugated derivatives appear in the equation we use the F 
and r operators: 



= -1 = -{C,^g )g 

Vt = c^ Vt = c\Ct9-')g 

and obtain the T and F operators with the following components: 

= - Cx r^. = - C X {Cxg"^)g 
f . = C* = Xt {Ct9-')9 

The above operators give in turn the currents: 

Jl = Sx't.x Jl = Sx'hx 

which are conserved according to the Proposition 3.2: 



(242) 
(243) 



(244) 
(245) 

(246) 
(247) 



provided the function x solves the equation ( p38| ), x' the conjugated equation: 



X' 



a" {c.9-')g+ " -c' a* {Ct9-')9 - a^ * 



(248) 
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and 5 is the symmetry operator for the conjugated equation. 

Following the Corollary 3.3 we can reformulate the currents ( p46| ) in order 

to obtain the conservation law including only the d derivatives: 

fj = 4 - C ( Jf ) Jf = Jt - Q{Jt) (249) 

The current J' obeys the conservation law: 

+ d'J'^ = (250) 

Now we apply the (^Tf ) version of the conserved current so as to obtain the 



conserved charges: 

Q' = j dxJf = jdx (^Sx'TtX - CiSx'ttx)) (251) 
where the integral J dx is defined as in (^). 



5 Final remarks 

We have discussed the method of derivation of the conservation laws for a 
class of equations with variable coefficients. It can be applied to models built 
using supersymmetric, discrete, mixed discrete and noncommutative (quan- 
tum Minkowski or braided) differential calculus. 

The conserved charges were constructed explicitly for considered supersym- 
metric and discrete equations. The problem of derivation of such charges 
for noncommutative models is open. The integration over noncommutative 
spaces namely braided linear space (including q-Minkowski) is well devel- 
oped [0, |, D, |14| but subintegrals and commutation rules for subintegrals 
and derivatives need further study. In particular the generalized version of 
the property ( p^ ) 

d' [ = [ d' 

J sub J sub 

must be derived. As we have shown for models on quantum planes [|14| it is 
crucial in the construction of the conserved charges. 
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6 Appendix A 



Proof of the Proposition 2.1: 

Let us recall the notation for monomials of derivatives: 



[pi...pk] := d''\..d'"^ \vx...v^\ := d -d (252) 

Due to the modification of the Leibnitz rule Leibnitz rule we are to consider 
the solution of the operator equation for the F operator in the form of the 
polynomial of order — 1: 

1=1 k=0 

where the coefficients depend on the coordinates x. 

The condition ( ^5]) from the section 2 applied to the above polynomial yields 
the equations for coefficients 

Af-l I 

1=1 k=0 II 

N-1 I 

+ 5Z t^l' H ^^+1' 
1=1 k=0 IJ, 

N-l I 

+ I] t^l' •••'^'^•] 5Z(^^^^'m...M;)l^fe+l' - , 
^=1 fc=0 M 

II II II 

= A(a) - K(d) 

The procedure is analogous to the one used in the proof of the Proposition 
4.1 from [0 which describes the derivation of conservation law for the equa- 
tion with constant coefficients or fulfilling the strong condition (see (23) and 
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(24) from [0). Now we have decided to change the form of the conjugated 
equation. This results in the weaker restrictions for coefficients of the equa- 
tion (|1]) and changes the set of equations for functions a^^j...^; defining the 
operator F for the following one: 



+ = A 



= 

"A, 



fi "'afii...fj.i 
k+1 



A, 



a/xi...^tfeAt/ifc+2---A»! 



"'lifil...fj.l I ^ "■afj.i...ij.iii 



with / = 1, - 1 



k 



0,...,l 



1. 



afii...fiktJ.tJ.k+i---fJ-i 
Oil (-- "1 A 



(255) 
(256) 
(257) 
(258) 
(259) 
(260) 



We see that the equations ( p55| , |256| , |260|) yielded by the coefficients of the 
conjugated equation differ from the case studied earlier in [|I2], |13[. We be- 
gin to solve this set of equations by deriving the coefficients a°^^ from 
equations ( |256| , p57| , p58| ). Namely for / = — 1 we have: 



S/i"'Q/ii.../i]v-l Wl---/iAr-i 

Applying the inverse operator we obtain: 



A,, 



(261) 



,0 



"'fifii...fiM-i ~ ^fj' A.Q,^^...^^_j (262) 
We insert this solution into ( P58|) for / = A^ — 2 and solve the next equation : 



an 



ti^afii...fj.N-2 



a/i/il.../i]V_2 



A 



/x^l.../i]V_2 



(263) 



By assumption (|]) after using ( operator and (^T]) we derive a' 





flfll...flM-2 



as: 



1° 

^J,fll...flN-2 



a^i...^iN-2 



(264) 



Passing to the next equation from the subset (p58D and solving them in the 
similar way we obtain the unique solution for coefficients a^: 



>l.../i; 



C "A 



afi2...lJ.i 



1 = 1, ...,N 



(265) 



This solution for initial coefficients allows us to evaluate the remaining ones 
using ( P59|J260|) , namely we obtain the coefficients after writing the subset 
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for k = and solving it the way we solved the subset ( |256| , |257| , P^5^ ) for 
a^. The result is unique and looks as follows: 



a 



1 



C "^AaiaM2.../., / = 1, iV - 1 (266) 



The same method applied to subsets of ( |259| , |260| ) for A; = 1, A^ — 2 produces 
the unique solution of the set of equations for coefficients in the form: 

~ C/xfc '"•••C/xi ^ Aai...afc«/ifc+i.../ii / = 1,...,A'^ — 1 (267) 

The derivation of the explicit formulas for unique solution of the coefficients 
of the operator concludes the proof of Proposition 2.1. 
Let us notice once more that in the derivation of coefficients for the 
operator the crucial factors were the properties of coefficients of the equation 
([T| - ^ which enabled us to solve the equation (P3| ) in the explicit form. 
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